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'^ ' Abstract 

- In this letter we consider models with N U{1) gauge fields A" together with 

rS I N Kalb-Raniond fields -B" in the large N limit. These models can be solved 

Cd ' explicitely and exhibit confinement for a large class of bare actions. The 

confining phase is characterized by an approximate "low energy" vector gauge 
symmetry under which the Kalb-Ramond fields -B" transform. A duality 
transformation shows that confinement is associated with magnetic monopoles 
condensation. 
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1 Introduction 

In various approaches towards a description of the confining phase of gauge 
theories - on the lattice, through duahty with the Higgs phase, ANO-strings or the 
"confining string" - the introduction of a Kalb-Ramond fields B^^^, [jl[| as effective 
variables has proved to be very useful [2 - 10]. In the case of pure Yang Mills 
theories they can be introduced as auxiliary fields for the abelian components of 
the field strength [8,9] in the maximal abelian gauge [^ and, after performing the 
path integral over the " non- diagonal" gauge fields (associated with non-diagonal 
generators), one is left with an effective action involving abelian gauge fields and 
Kalb-Ramond fields only. 

Common to all these approaches is the idea that Kalb-Ramond fields are effective 
variables only at low energy, i.e. in the infra-red regime. It is thus natural to provide 
models with Kalb-Ramond fields as effective variables with an ultra-violet cutoff, 
but to allow for irrelevant operators in the corresponding bare action. On the other 
hand it is sufficient to restrict oneself to abelian gauge theories (or the abelian 
subsector of non-abelian theories): Monopole condensation, which is believed to be 
the mechanism behind confinement [12, 13], is a purely abelian phenomenon. In 
the context of non-abelian gauge theories the relevance of the abelian sub-sector is 
conventially referred to as " abelian dominance" . 

The purpose of the present paper is the study of models with N abelian gauge 
fields A^ and A^ Kalb-Ramond fields 5" in the large N limit. As we will see, they 
can be solved using standard functional methods employed for large N field theories. 
The emergence of a confining phase can be seen explicitely, and the appearance of 
a "low energy" vector gauge symmetry allows for a duality transformation showing 
that magnetic monopoles have condensed. In this letter we present the essential 
results of this approach, leaving many details (as the precise relation to Yang-Mills 
theories) to a subsequent publication [|l4 . 



The A^ - B^^y — model 



The starting point is the partition function for a model with the above field 
content. Adding sources for A" and B^^^, n = 1 . . . N, and a covariant gauge fixing 
term the partition function reads 



Due to the N f/(l) gauge symmetries Sharei,^, B) can only depend on F^y, hence 
we can write Sbare{F, B). Next, in order to allow for a large N expansion, we assume 
that Sbare{F,B) dcpcuds on 0{N) invariants (singlets) only. The aim would be to 
allow for a dependence of Share{F, B) on 0{N) singlets as general as possible. Here 
we confine ourselves to the following ansatz: First we introduce three Lorentz scalar 
0{N) singlet operators 

N 

n=l 

N 

02{x) = j:F;:,ix)B;,ix) , 

n=l 

N 



o,{x) = Y.b;,{x)b;^{x) . (2.2) 

n=l 

Then we take Sbare{F, B) of the form 

Sbare{F,B) = jd'x \Cbare{0,) + ^ (5,5;,)' + | (^M^?;.)'} (2-3) 

where 

B^^u = 2 ^A'^P'^ 5p<7 • (2-4) 

We allow Cbare in ( |2.3| ) to contain arbitrary derivatives acting on the operators Oi. 

This is still not the most general form of Sbare{F, B); one can certainly construct 

infinitely many more 0{N) singlet operators which contain open Lorentz indices 

and/or "internal" derivatives as the second and third terms in (2J). It can be 



argued ||1J| , however, that these do not modify the essential features of the results 
obtained below. 

In order to solve the model in the large N limit we have to make assumptions on 
the N dependence of the parameters in Sbare- These assumptions can be summarized 
by rewriting (|2.3|) as 

SbareiF, B) = J d'x ^N Cbare Q) + ^ (^M^;.)' + f (^.^M.)'} ^^.5) 

where now the coefficients of Cbare are independent of A^. 



3 The large N solution 



The most convenient formalism for the treatment of field theories in the large 
N limit is the introduction of auxiliary fields for composite 0{N) singlet operators 
||15|| . In the present case we introduce one auxiliary field (pi for each of the bilinear 
0{N) singlet operators Oi in eq. (|2.2|). This amounts to re- write the term involving 



NCbare (C'j/A^) in the exponent of ( ^.11 ) as 

^-Njd'^xCtarei^) ^ ]_ f -J^J). ^-NGtarei'P.yJd'^X^.O, . (3.1) 

In the large A^ limit the path integral on the right-hand side of (|3.1|) can be replaced 
by its stationary point, and the relation between Gbare and Cbare becomes 

TV J d^X Cbare (^) = NGbare{ct>i) + / d^X 0, O, . (3.2) 

Equation (|3.2|) allows, in principle, to construct Gbare{<Pi) from Cbare, although here 
we allow Cbare to be an arbitrary functional (including derivatives) of Oi. Next we 
insert eq. (|3.1|) into (|2.1|), which becomes 



(3.3) 

The VA VB path integrals have become Gaussian in (|3.3|) . In order to express the 
result in compact form we introduce the notation </?" = |A",i?" >, i.e. the indices 
r attached to the fields </?"■ denote both the different fields A", B"^ and the different 
Lorentz indices. Correspondingly we introduce notation J" for {Ja,^i, Jb,^iu}- The 
result of the Gaussian integration over VA T>B can now be written as 



g-iy(J) ^ J_ /"-p ,. ^-NGtare{<^>^)-NAG{(|>^) + ^ J d'ixld'ix2{J?(x^)P^''{x^,X2,<^>^)J^{x2)} (^^ ,^\ 



with 



AG'(0.) = -^ Tr log (P^^(xi,X2,0.)) . (3.5) 



The propagators P^'* of the A^, B"^^ - system are proportional to 5n,m with 
n,m = 1 . . .N and we took care of the resulting contribution from the trace in ( p.5|) 



by the explicit factor A^ multiplying AG in ( plil) . The propagators -P^^, Papa ^^d 



^l,pa 



Pj^i^pu depend on the terms 0iCj, |(9-B)^, |(9i?)^ and -^{dAf in the exponent of 
( p.3|) . Simple explicit expressions can be obtained only for constant fields (pi] in this 
case one finds for AG (in the Landau gauge a — > 0) 



AG((/., 



d X 



d!^p 
(2^ 



log(0iOrp2 + 4(/.i(/)3 - (g) + log (V + 



(3.6) 



The d'^p integral in (p.6|) has to be performed with an UV cutoff A^. The result 
simplifies considerably if one introduces 



H(Pz 



acpi 



(3.7) 



Up to field independent terms one then obtains 



AG(0,) 



327r2 



/ d^x [(A^ - S2) log (a0i(A2 + S)) + T? log ((T0iS) + A^S 



+ fA4_160i^,_^,2, , .^A , 1602 ^_^^^^ , ^^2<^3 



/l2 



log (A^/i + 403) + -^ log (403) + 4A 



h 



(3.8) 



This expression for AG(0i) has to be inserted into ( |3.4| ) and, in the large N limit, 
the 'Dcpi path integral is again dominated by its stationary point (s). Hence W{J) 
becomes 



W{J) = NG{(Pi 



^Jd^Xid'x2{j:{Xi) P'^'{X,,X2A^)J:{X2)} (3.9) 



where 



and 



G'(0.) = G,,„(0^) + AG(0,) 
i{J) satisfy the three equations (recall i = 1, 2, 3) 



(3.10) 



■(ArG(0,) -^Jd'x,d^X2 Jr(xi) P'-'{x,,X2,(P.)J:{x2)) 



. (3.11) 



MJ) 



The model is thus solved, for given Gbare{(pi), up to the technical problem of finding 
the stationary points 4>i{J). 



Next we wish to show that the particular configuration where 



■>W3 



2 = . (3.12) 

(or S = 0) is a "natural" solution of the three stationary point equations ( |3.11|) , i.e. 
a solution which requires no fine tuning of the parameters in Gbare{4'i)- Below we 
will see that the phase where ( p.l2|) holds is the confining phase of the model. 

In order to "see" the solution ( p.l2| ) of the eqs. ( |3.11| ) it is necessary to regularize 
the singularity of the derivatives of AG(0j) w.r.t. the fields at S = 0. The origin 
of the non-analytic behaviour of AG{4>i) at S = is the infrared behaviour of the 
propagators P^''{(f)i) which, in momentum space, behave like P^^{q'^,(j)i) ~ g~^ for 
S = 0. In order to regularize these infrared singularities we perform the d^p integral 
in ( |3.6|) also with an infrared cutoff k"^. For J" = (hence we write 0? instead of 



) the three stationary point equations (|3 . 1 1|) can then be brought into the form 



6G. 



bare 



+ 




Hi 

p^ + h^ 



32tt^c 



(a" - k' 



^M^)^^'-'' 



J J Jo 







(3.13) 



The solutions 0° of ( |3.13|) , and hence S°, depend on the infrared cutoff k'^. Let us 
now assume that 



6G 



bare 



sj: 



S0=0 



167r2 



A^ >0 



(3.14) 



Then the last of the stationary point equations (|3.13| ) implies S*^(A;^) < for 
A^ > /c^ > 0, and Tj^{k'^) behaves as follows for fc^ -^ 0: 



S°(A:^ 



0_ 



16n' 



-S°log 



S0 + A2 



5G, 



bare 



sj: 



S0=0 



167r2 



A^ . (3.15) 



Hence, under the condition ( p.l4| ), we obtain S(0) =0 naturally. Note that this 
stationary point would not have been observed if one puts k^ = from the start. 

What is the meaning of the condition ( p.l4| ) or, better, S°(A;2) < 0? To this 
end we push the infrared cutoff k"^ upwards until it reaches the UV cutoff A^, and 



investigate the consequence of E°(A^) < . Given the definition (|3.7| ) for S, and 
for 0?(A^) > 0, this latter condition reads 

40?(A2)0°(A2) - {^l{A')f < . (3.16) 

For k"^ = h? the contribution AG{(j)i) to G{(j)i) in (|3.10| ) vanishes, and the configu- 
rations 0°(A^) are the stationary points of Gbare{<Pi)- From the Legendre transfor- 
mation (|3.2| ) and the definitions ( p.2|) of the operators Oi it is now straightforward 
to see that the inequahty ( p.l6|) imphes (with ys" as below (|3.3| )) 



j-^ , / " *~'bare 

Vet 



< . (3.17) 

F=B=0 



( p.l7|) corresponds to a ChareiF", B) which is non-convex at the origin of field space. 
All this is very similar to the case of a non-convex bare scalar potential: The 
effective scalar potential has to be semi-convex, and if the bare scalar potential is 
sufficiently non-convex the effective potential becomes flat in its "inner" region. As 
in the present case the observation of this phenomenon requires the introduction 
of an "artiflcial" infrared cutoff fc^, and a careful discussion of the limit fc^ — > 
|T^. We emphasize that consequently the emergence of the "conflning" phase ( |3.12| ) 



(see below) is not a particular feature of the large N limit; the advantage of the 
large N limit is only to allow for an explicit study of this phenomen (for given 
parametrizations oi Cbare{F,B) or Gbare{4>i))- 

4 Properties of the confining phase 

In the following we assume that the necessary inequality on the parameters of 
^bare{F,B) (a Sufficiently negative curvature at the origin of field space) for the 
reach of the confining phase ( |3.12| ) is satisfied. In order to discuss its properties it 
is more convenient to switch from W{J) in (|3.9| ) to the effective action T{A, B) via 
a Legendre transform with respect to the sources J: 

r{A, B) = W{J) + J d'x {JI^AI + J5,,,5;,) . (4.1) 

As discussed in detail in |T^ one obtains 

r(A, B) = NG{$,) + J d'x Uo, + I [d.B^^Y + I {d,B;,f\ (4.2) 



with G{(j)i) as in ( p.lO| ), and the stationary point equations (|3.11|) for 0, can be 
written as 



First we note that because of the relation ( p.l2|) the expression (piOi in ( ^l2|) 
becomes 



o^ = J:[ V 01 C + V03 b;^ 



(4.4) 



Consequently it is invariant under the following gauge symmetry involving vector- 
like gauge parameters A^ |jl[]: 



SAlix) = K{x) , 5F-{x) = d,K{x) - dAlix) = A" (a:) 



6B;{x) 



\ 



03 



(4.5) 



In addition one finds that the term 



{d^B^,.? m 



21) is also invariant under 



( [4. 51) thanks to a Bianchi identity, provided the configurations (pi are constant in x. 
(Note that, from eq. ([4.3|) with F as in (^4.2| ), constant configurations (pi result from 



constant configurations Of, however, from their definition (|2.2| ), constant Oi do not 



necessarily imply constant configurations F" and -B" .) The last term ~ [d^j^B^^Y 
in ([4.2|) behaves as a gauge fixing term of the symmetry ( [4.5|) , and its presence 
insures the existence of the propagators. 



It is to be expected that the symmetry ( [4.5|) is broken by higher derivative 
terms (beyond the gauge fixing term): The bare action Si,are (|2.5| ) of the model 
does certainly not exhibit the symmetry (14.5|) , and the Green functions at large 
non-exceptional Euclidean momenta with p^ -^ A^ are generated by Share- This fact 
is realized by the dependence of the effective action on higher derivative terms. The 



symmetry ( |4.5| ) is thus a pure "low energy" phenomenon. The implication of the 
gauge symmetry ( ^I5| ) on modes of the U{1) gauge fields A^ which correspond to 
constant configurations Oi is that they can be "gauged away" and "eaten" by the 
(massive or even infinitely massive) Kalb-Ramond fields -B" , in some analogy to 
the ordinary Higgs effect |]I[. 



In addition the symmetry ( [4.5|) allows for a duality transformation: A priori the 
dual of a U{1) gauge field A^ (in d = 4) is again a. U{1) gauge field C^ (whose field 



strength tensor will be denoted by -F^;^)? ^^id the dual of a Kalb-Ramond field 5" 
is a (pseudo-) scalar </?". In the present case the duality transformations mix the 
fields and read 



^F;;r = ^if;, + ^fM^B;, 



9m^" + C 



h 

2\ 



t O^B:^ (4.6) 



where the tildes on B^^ and F" have been defined in (|2.4|) . The corresponding dual 
action reads (as obtained from ([4.2|) without NG{(j)i) and without the "gauge fixing" 
term) 

r,,..Kc^,¥') = /rf'^|^^;;r^;;r + ^|^(5.¥'" + c^;)'} • (4.7) 

Note that, due to the implicit dependence of 0j on F and 5, these duality 



transformation are non-linear. Actually one finds |TJ] that only half of the equations 
of motion and Bianchi identities are exactly interchanged through ( [4.6|) and ( [4.7|) , 
whereas the other half holds again only for constant configurations 0j and hence (9,. 
Thus duality is realized at the non-linear level again only in the corresponding " low 
energy" regime. 

The physical interpretation of the dual action ( ^4.7| ) is obviously the one of an 
abelian f/(l)^ Higgs model in the spontaneously broken phase where yj" repre- 
sent the Goldstone bosons, and where the gauge fields C^ have acquired a mass 
2(03//i)^/^. Since this represents the "low energy effective action" of a theory in 
which the "dual" electric charge has condensed in the vacuum, the original action 
( [4.2|) with (|4.4| ) corresponds to the situation where the "magnetic" charge has con- 
densed in the vacuum. 

Let us turn to the response of the model in the confining phase with respect 
to external sources. The expression for W{J) has been given in eq. (|3.9|) in the 
preceeding section, and first we concentrate ourselves on the term quadratic in the 
sources J". Let us start with a source Ja,^{x) for the fields A^ only, which is of the 
form of a Wilson loop: 



JaA^) = i9A J dx'^ 6 {x - x') . (4. 



The term quadratic in J in ( p.9|) then becomes 



j dxi,^ j dx2,u P^,u i^i - ^2) (4.9) 



where P^^ has to be obtained from the action ( ^.21 ) with ( |4.4| ): 



P^,^(') = 7T^ 1 V I M^ - ^ l°g 1^1 + ^°^«t. 



I67r^0i t '^ Vl^r ^ 



-^9^9jlog|;2|-^^(log|^|+const.')j| . (4.10) 

The (actually divergent) constants in (|4.1CI|) disappear in the expression (|4.9|) . In 
the limit where the (minimal) surface S enclosed by the loop C in ( [4.8[) becomes 
very large one finds that the expression ( [4.9|) is proportional to S, thus one obtains 
the area law for the expectation value of the Wilson loop. 

However, at first sight an inconsistency arises due to the long-range behaviour 
of the propagator P^^^{z) in ( [4.10| ): Let us imagine that space-time is filled with 
"virtual" Wilson loops (originating, e.g., from vacuum bubbles of virtual quark- 
antiquark pairs), and let us compute the corresponding contribution to the action 
due to the interactions among different "virtual" Wilson loops. Even if one assumes 
that these "virtual" Wilson loops are arbitraryly tiny in size, localized in space- 
time and if one averages over their orientation in space-time, the contribution to 
the action induced by the long-range behaviour of the propagator P^^^{z) in ( ^.10| ) 
diverges (logarithmically) in the infinite volume limit. 

This infinity can be avoided, however, once one realizes that all components 
{r,s} = {A,B} of the propagators P^'^lz) in ( |3.9| ) have a "bad" long range be- 
haviour. The precise expressions for all propagators, as obtained from the action 
( [4.2|) , will be given in |T^. One finds that all terms in the propagators which de- 



crease not sufficiently fast at infinity in order to avoid the above infrared divergence 
(which originate from g~^-terms in momentum space) cancel in the sum over r and 
s in J^P^^J^ in (|3^ ) if the sources J^^^ and Jb,^u satisfy 

OsJl^ix) = 2^/$^^,Jl,^{x). (4.11) 

This result can also be phrased as follows: In the presence of an arbitrary back- 
ground of "virtual" Wilson loops it costs an infinite amount of action (or energy) to 
"switch on" sources J^^^ and/or J^^^u "which are not related as in ( [4.11D , due to the 

10 



interactions induced between the sources and the background of "virtual" Wilson 
loops induced by the long-range terms in the propagators. 

If the source J^^ is of the form of a Wilson loop ( [4 .81) one finds that ([4.11|) 
implies that the source J^^^^ is of the form of a " Wilson surface" , 

Jb,,Ax) = WB I d'a,,{z) 6\x - z) , (4.12) 

where the surface S is bounded by the loop C in ( [4.8|) (but otherwise arbitrary) and 
where gs satisfies y 0i gs = y 03 5'a/2. 

It is straightforward to see that the condition ([4.11|) on the sources is equivalent 



to the condition that the couplings J^^^^J] + JB,fiuB^u o^ ^^e fields to the sources 
respect the "low energy" gauge symmetries ( [4.5| ). In the case of conventional gauge 
symmetries these conditions can (and have to) be imposed by hand in order to 
ensure renormalizability and unitarity of the theory. In the present model, on the 
one hand, they cannot be imposed from the beginning, since the associated (low 
energy) gauge symmetries appear only at the level of the effective action once the 
equations of motion of the fields 0j are satisfied. Although renormalizability is not 
an issue here, since we consider an effective low energy theory with a fixed UV cutoff, 
it is interesting to see that the corresponding condition on the sources is generated 
dynamically in the sense that its violation costs infinite action. 

Clearly we now have to reconsider the expectation value of the Wilson loop, 
which consists now of a source (|4.8| ) for A" and a source ( [4.12| ) for i?" inserted into 



the term quadratic in J in ( |3.9|) (higher orders in J will be discussed below). Using 
all propagators P^^ from |T^ this term becomes 





s Js le-K^cpi V h 

X (T^,^u,pAd) - ^T2,^,,pJ \z,-Z2\ ^' I ^'-"^ " ^''^/^ I ^^-^^^ 

with 



Ti,^u,pa{d) = 5^pdud„ - Spadudp - S^p dpd„ + 5y„ dpdp , 

-'■2,^,1', pa "/^p "j/CT ('p.a ^up 5 V / 



11 



and where Ki is a Bessel function. In the hmit where the surface S becomes large 
the expression (|4.13| ) behaves as 



S.^^r^]''' TdzK^zM . (4.15) 

7i(Pi\hJ Jo V / 

Hence it imphes the area law in spite of the cancellations of the long range contri- 
butions of the propagators. (Since we had omitted the UV cutoff in the space-time 
propagators the expression ( |4.15| ) is seemingly UV divergent). 

The preceeding results, based on a treatment of the term quadratic in J in (^.91), 
have obvious interpretations in the context of the stochastic vacuum model [|T^ for 
Yang-Mills theories: There, in the Gaussian approximation, the expectation value 
of the Wilson loop is given by the expectation value of the field strength correlator, 
which plays the same role as the term quadratic in J in (|3.9| ). At first sight an 
ambiguity appears: A priori it is not clear, whether the Yang-Mills Wilson loop 
reappears in our "effective low energy model" (after integrating out the off-diagonal 
gluons in the MAG, see fl^) in the form of a source Ja, Jb or, most likely, as a 
combination of both (or even in the form of additional terms in Share)- The condition 
( [4.11| ) fixes this ambiguity. 



If we identify naively the term quadratic in J in ( p.9|) with the expectation value 
of the field strength correlator, the area law obtained from ( [4.9| ) above corresponds to 
a function Di(x^) in the standard decomposition of the field strength correlator fT?]) 
which decreases only as \x\~'^ for large |a:|. Such a behaviour is strongly disfavoured 



by lattice measurements |[TH| of the Yang-Mills field strength correlator. On the 



other hand the results ( ^.13| ) and ( [4.15| ), which are the consequence of the condition 



( [4.11[ ), agree well with with the lattice measurements of the Yang-Mills field strength 



correlator [1131 and are in fact identical to the results for this correlator obtained in 



various models 10, IS, 



We recall, however, that up to now we have only discussed the term quadratic in 
J in (|3.9|) . Let us note that this term coincides (up to an irrelevant constant) with 
W{J) to 0{J'^), provided that we replace 0j by 0^ in P^^{(f)i): From the stationary 
point equations (|3.11|) for 0j we have 4>i{J) = 0° + (9(J^), and - since the 0° are 
stationary points of NG{(j)i) - we thus have NG{(j)i{J)) = NG{$'^i) + 0{J^). Beyond 
an expansion in powers of J the stationary point equations (p.ll| ) are, however, 
cumbersome to solve since they involve the full propagators P'"'^(0j). 

To this end it is wiser to start with the effective action T{A, B) as given in eq. 

12 



^. One should to solve the combined equations of motion for A" 5" in the 
presence of sources J^^^ and Jb^^^u together with the equations ( [4.3|) for 0j (in some 



analogy with the approach in |T9| based on the dual abelian Higgs model). The 
corresponding solutions have to be inserted into T{A, B) in ([4.2|) , and then one has 
to "undo" the Legendre transformation ([4 .11) in order to obtain W{J). This last 
step is actually trivial since T{A,B) is quadratic in A, B: It suffices to change the 
sign of the second term in the expression ( [4 .21 ) . Then one can study the dependence 
oiW{J) on J in its full beauty. 



As discussed in somewhat more detail in |]n[ the effect of such a more complete 
calculation can be estimated in the simple case where the sources J are non- vanishing 
only inside a finite volume V (to be identified with a "Wilson surface" of finite width 
AL), and where derivative terms in G{(f)i) are neglected: Then the solutions (j)i{x) of 
the stationary point equations differ from 0° inside V , but coincide with 0° outside 
V . Hence one obtains an additional contribution to the action proportional to V , 
or a contribution to the energy of the configuration proportional to the diameter 
of the Wilson surface (at fixed time t). This picture supports "fiux-tube" models 
for the origin of the string tension. Its details (such as the shape of the fields (pi 
perpendicular to the string axis) depend, however, on the precise form of Gbare{(pi) 
and hence of Cbare{Oi). (At this point it is useful to note that the auxiliary fields 
(pi, as introduced in ( |3.1| ), parametrize vevs of the bilinear operators Oi in (|2.2| ).) 

We close this section with a comment on the physical spectrum of the model. 
In momentum space the propagators P^''(g^) for the A^-B^^ - system have, in the 
confining phase, q~^ singularities for g^ ^ as well as poles at q^ = —A(f)3/h. The 
q~^ singularities for g^ ^ do not correspond to asymptotic physical states. The 
poles at g^ = —icps/h, on the other hand, would disappear if we would replace the 
constant h in the ansatz ( |2.3| ) for Share by a function /i(g^) such that /i(g^) vanishes 
sufficiently rapidly for large |g^|. This would not modify any of our essential results, 
but would be motivated by the idea that the Kalb-Ramond fields 5" have originally 
be introduced into a more " microscopic" theory (as a Yang-Mills theory) as auxiliary 
fields (for, e.g., the abelian field strengths F^^, P, ^ Q), and their kinetic terms 
are thus loop-effects of modes which have been integrated out (as the off-diagonal 
gluons). Then none of the degrees of freedom of the A^-B^^ - system would appear 
as asymptotic states, consistent with the absence of {N^ — 1) - plets in a SU{Nc) 
Yang-Mills theory. 

We have also searched for "bound states", i.e. poles in the propagators of the 0j 
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fields. Such poles are not present, essentially because bosonic loop contributions to 
kinetic terms of auxiliary 0{N) singlet fields differ in sign with respect to fermionic 



loop contributions (which do generate propagating bound states ||T5l). One may 
be deceived by the absence of "glueballs" , if the present model is interpreted as an 
effective low energy theory for SU{Nc) Yang-Mills theory. However, we recall that 
the present model would only describe the abelian subsector of SU{Nc) Yang-Mills 
theory in the MAG p, ^, Q) and that we finally have to add the off-diagonal gluons 
as well. Our theory induces confining interactions among all fields which couple A^ 
and B^j^, hence the off-diagonal gluons will necessarily form bound states which will 
correspond to the desired glueballs. 

Thus, if we use functions /i(q'^) with the above-mentioned properties in the ansatz 
]3|) for Share, the model has no asymptotic states at all. Its only "meaning" is then 



to react to external sources, and to confine them as discussed before. 

5 Discussion and Outlook 

We have studied a class of four-dimensional f/(l) gauge theories including Kalb- 
Ramond fields, which exhibit confinement and allow nevertheless - in the large 
N limit - for controllable computations in the infrared regime. Some features of 
the confining phase correspond quite to our expectations, notably the possiblity 
to perform a duality transformation of the low energy part of the effective action 
and thus to interpret confinement as monopole condensation. A technically related 
phenomenon is the appearance of a low energy vector gauge symmetry, which allows 
to "gauge away" the low momentum modes of the abelian gauge fields A^. 

An interesting feature is the origin of the relation among the parameters of 
the effective action which generates the above symmetry: The bare action, as a 
functional of A" and 5^,^, has to be (sufficiently) non-convex, which renders the 
effective action "fiat" in some region around the origin in field space. This "flatness" 
corresponds to the above symmetry with all its consequences. This phenomenon is 
evidently independent from the large N limit employed here. The auxiliary scalar 
fields (pi parametrize various bilinear condensates which do not, however, break any 
internal symmetry. 

We have argued - but not shown in detail (to this end see |]I4[) - that the present 
class of models can be obtained from SU{Nc) Yang-Mills theories in the maximal 
abelian gauge after integrating out the off-diagonal gluons. This application clarifies 
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why it is sensible to consider these models as effective low energy theories equipped 
with an UV cutoff, but with a bare action including non-renormalizable interac- 
tions. Our ansatz (|2.3| ) for the bare action is already quite general and exhibits the 
most interesting phenomena, but it could easily be generalized by including further 
bilinear operators with "internal" derivatives and/or external Lorentz indices. It 



can be argued |IJ] that further bilinear operators with "internal" derivatives do not 
affect the low energy limit of the model (and just modify somewhat the relation 
between the bare and effective actions), whereas operators with external Lorentz 
indices will have no vevs (but can affect the response of the model with respect 
to external sources). More detailed investigations in this direction would be quite 
st r aight f or war d . 

A grain of salt constitutes the fact that the large N limit in the present class of 
models does not coincide with the large Nc limit of SU{Nc) Yang- Mills theories fl^ . 



First, the powers of N^ in Share would not correspond to the powers of N required 
in (|2.5|) , and second Share would not necessarily depend only on 0{N) singlets since 
this 0{N) is not a sub-group of SU{Nc) (with Nc = N + 1) . A discrete reflection 
symmetry, under which A^ and -B" (or the abelian components of F^iy) change sign 
is, however, a symmetry of SU{Nc) Yang-Mills theories, which justifies at least the 
introduction of bilinear composite operators. We recall again, on the other hand, 
that the essential features of the confining phase do not rely on the large N limit. 

We have seen that the "physics" of the confining phase is by no means unique. 
On the one hand confinement can always be interpreted as monopole condensation 
(and hence the vacuum as a dual superconductor), but many features like the most 
important contributions to the string tension, (non-local) vacuum correlators and 
the shape of the vacuum energy distribution perpendicular to a flux tube depend 
on the non- universal properties of the model incoded in Share- Hence, if we wish to 
learn more about the way confinement is realized in Yang-Mills theories we have to 
find ways to learn more about the bare action, and eventually to handle the present 
class of models beyond the large N limit. Nevertheless the solvable version of the 
present models can certainly play the role of a useful laboratory for the study of the 
properties of a confining phase in the future. 
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